r e b l i q u e 



FRANCA 



10/50 13 ?Q JJsPI 

Rec'd PGT/PT© i 4 JUL 2004 



IN8TITUT 
NATIONAL DB 
LA PROPRIETE 
INDUSTRIELLB 



^€I/F^ 0 3 / 0 0 0 9 9 




BREVET D'INVENTION 



CERTIFICAT D'UTILITE - CERTIFICAT D'ADDITION 



COPIE OFFICIELLE 



Le Directeur general de l'lnstitut national de la prpprlete 
industrielle certifie que le document ci-annexe est la copie 
certifiee conforme d'une demande de titre de propriete 
industrielle deposee a I'lnstitut. 



Fait a Paris, le 



2 1 JAN. 2003 



DOCUMENT DE PRIORITE 

' PRESENTE OU TRANSMIS 
| CONFORMEMENT A LA 
i REGLE 17.1. a) OUb) 



Pour le Directeur general de I'lnstitut 
national de ta propriety industrielle 
Le Chef du DSpartement des brevets 




Martine PLANCHE 



BEST AVAILABLE COPY 



SIEGE 

I N S T I T U T 26 bis, nio do Saint Petersbours 
NATIONAL DE 75800 PARIS ccdox 08 
national ut Telephone : 33 (1) 53 04 53 04 

LA PROPRIETE Tdldcbplo : 33 (1) 42 93 59 30 

INDUSTRIELLE vww.lnpLfr 



ETABLISSEMENT PUBUC NATIONAL CREE PAR LA LO! N* S1.4AA mi is AV »tl iori 



* 



■ 1HST1TUT 

LAFffomiKTB 
I rtOUflTRI ELLQ 



BREVET P'IMMENTiON 
CERTIFICAT C^PLITE 

Code de la propriete intellectuelle - Uvre VI 



N° 11354*01 



26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 33 (1) 53 04 53 04 Tetecopie : 33 (1) 42 94 86 54 



REQUETE EN DEUVRANCE 
page 1/2 



Get imprime est a remplir lisiblement a Tencre noire db wow/3co3oi 


REMISE DES PIECES 
DATE 

UEU 14 JAM 2002 
75 INPI PARIS B 

"° O'ENREGISTREMENT 0900^7 
NATIONAL ATTRIBUE PAR UlNPI %f£Z.**P*U>^^ ff 

DATE DE DEPOT ATTRIBUEE 

i % jAN.-am 


g§ N0M ET ADRESSE OU DEMANDEUR OU DU MANDATAIRE 

A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 
□ ■ 

n * ■ 


Vos references pour ce dossier 

(factdtatifi 


Confirmation d'un d6p6t par telecopie 


□ N° attribuS par I'INPI a la tetecopie 



Demande de brevet 



Demande divisionnalre 


□ 




Date 1 1 1 : 1 1 1 1 .1 
Date III 1 i i _J J 

Date 1 i 1 i 1 i | i | 


J Demande de brevet in itiale 


N° 




ou demande de certificat dutilite in itiale 

Transformation d'une demande de 
brevet eu ropeen Demande de brevet initiate 


N° 

□ 




|H TITRE DE L'INVENTION (200caractdres ou espaces maximum) 


oe 













DECLARATION DE PRIORITE 
OU REQUtTE DU BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANT&R1EURE FRANQAISE 



Pays ou organisation 
Date I | I » 1 i 



Pays ou organisation 
Date hlil. 



N° 



Pays ou organisation 

Date 1 i I i 1 i < ■ I 



□ S'il y a d'autres prio rttes, cochez la case et utttisez rimprimS «Suite» 



Nom ou denomination sociafe 



Prenoms 
Forme juridique 
N°SIREN 
Code APE-NAF 



j ? „ l J, ,«, , l „ L.L l J 



Adresse 



I i i i 



Rue 



4 ;Uuu^vJ:O p 4t6^^ic 



Code postal et ville 



Pays 
Nationality 

N° de telephone (facidtajjf) 
N° de til(§copie (facultaiij) 
Adresse electronique (facidtatif) 



Remplir imperativement la 2 6w page 




• 



BREVET D'BiMMfcFOlO^ 



CERTIFiCAT D'OMTE 

Code de la propriete intellectuelle - Livre VI 



N° 11354*02 



% rnopmfji 



REQUETE EN DELIVRANCE 

page 1/2 



Cet imprime est a remplir lisiblement a I'encre noire 



inDUSTRlELLC 

26 bis. rue de Saint Peterobourg 
75800 Pans Cedex 08 

Tcleohcme : 33 (1) 53 04 53 04 Tulccopie : 33 (!) 42 94 86 54 



DO Ml) V// 010501 



REMISE PES PIECES 
DATE 

LIEU 

N* D'ENREGISTREMENT 
NATIONAL ATTRIBUE PAR L'INPI 
DATE DE DEPOT ATTRIBUEE 
PAR UINPI 



A Reserve a IMNPl] 



A h- 1 ■ o z 

o%oo S3 ? 



Vos references pour ce dossier 

(facultaiif) 



239561 D19929 JR 6- 



| NOM ET ADRESSE DU DEMANDEUR OU DU MANDATA1RE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 

Cabinet REGIMBEAU 
20, rue de Chazelles 
75847 PARIS CEDEX 17 
FRANCE 



Confirmation d'un depot par telecopie 



Demande de brevet 



Demande de certificat d'utilite 



Demande divisionnaire 

Demamle de brevet initiate 
ou demande de certificat d'utilitt initiate 
Transformation d'une demande de 
brevet europeen Demande de brevet initiate 




T1TRE DE LMfWEfSTlOW (200 caracttres ou espaces maximum) 
PROCEDE DE FILTRAGE DE DONNEES, NOTAMMENT DE TYPE SISMIQUES, PAR ANALYSE KRIGEANTE. 



fU DECLARATION DE PR10RITE 
OU REQUETE DU BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANTEKIEURE FRAWQAISE 



Pays ou organisation 
Date | i 1 i I i i i I 
Pays ou organisation 
Date 1 i 1 ! I 
Pays ou organisation 
Date 1 i 1 i I i i i 1 
□ S'il y a d'autres prlorites, cochez la case et utilisez r imprime «Suite» 



Nom 

ou denomination sociale 



GQMPAGNIE GENERALE-DE- GEOPHYSIQUEr - 



Prenoms 

Forme juridique 
N° SIREN _ 
Code APE-NAF 



Rue 



Code postal et ville 



Domicile 
ou 

siege 

"Nationalise 

N° de tele pho ne (fa cultaiif) 

. Adresse electronique (facultaiif) 



Pays 



□ • Personne morale 



□ Personne physique 



^6^2Q22|H 



1 , run IgfiQ ¥Mf"T 21300 massy 



FRANCE 
"Franfaise 



N° de telecopie (facullatif) 



PI «'SI *' o nine A'ttn 




BREVET D'i 
CERTIF8CAT 



ENTI0N 
LITE 



5 1NCT1TUT 
KMIONAL OS 
LA PQOPaUTE 

mouoTKicttd 



REQUETE EN DELiVRANCE 
page 2/2 



REMISE OES PIECES 
DATE 



|R6serv6&riNPlh 



LIEU 



14 JAM 2002 
75 INPI PARIS B 

H° D'ENREGISTREMENT 
NATIONAL ATTRIBU£ PAR LT1NPI 



Vos references pour ce dossier : 

(Jaaittatif) 



DB 540 w 



Nom 
Prenom 

Cabinet ou Societe 

N °de pouvoir permanent et/ou 
de lien contractuel 

Rue 

Adresse 

Code postal et ville 
N° de telephone (facnftatij) 
N° de telecopie (JaailMiJ) 
Adresse electronique (facultatij) 



' <:, -f u*i w-.; ii^ 



Les inventeurs sont fes demandeurs 



□ Oui 

E Non Pans ce cas fournir ung designation d r inventeur(s) sgpgrjg 

mmi pour arte *d$rri$i?fie de %£y$t;(yr compns dhrtsJon et fr^nsfoHriatidri) 



, 



Etablissement immediat 
ou etablissement differe 



□ 



Paiement echelonne de la redevance 



Paiement en deux versements, uniquement pour les personnes physiques 

□ Oui 

□ Non 



REDUCTION DO TAUX 
DES REDEVANCES 



Uniquement pour 5es personnes physiques 

□ Requise pour la premiere fois pour cette invention (joindre tin avis de non-imposition) 

□ Requise anterieurement a ce depot (joindre une capie de la decision dad mission 
pour cette invention ou indiquersa reference): 





Si vous avez utilise I'imprime «Suite», 
indiquez ie nombre de pages jointes 






EE SIGNATURE DU DEftiANDEUR 

OU DU MANDATAIRE Ok/ //O? 
(Worn etqualitedu signature) 1 ^ vUL/ C 


VISA DE LA PREFECTURE 
OU DE L'INP! 



La loi n°78-17 du 6 janvier 1978 relative a Tinformatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de PINPI. 



# 

i 



PROCEDE DE FILTRAGE DE DONNEES. NOTAMMENT DE TYPE 
S1SMIQUES. PAR ANALYSE KRIGEANTE 

DOMAINE TECHNIQUE GENERAL - RAPPELS SUR L'ANALYSE KRIGEANTE 

5 

La presente invention est relative aii filtrage de donnees par 
analyse krigeante. 

Elle trouve avantageusement, mais non limitativement, application 
dans le domaine de la sismique. 
10 L'analyse krigeante permet la decomposition d'une fonction 

ateatoire a partir de sa fonction de covariance. 

Elle est notamment classiquement utilisee en geostatistique pour 
filtrer des donnees sismiques, notamment, mais non limitativement, pour 
caracteriser des reservoirs. 
15 L'analyse Krigeante repose notamment sur i'hypothese qu-un 

phenomene mesur6 localement au moyen d'un echantillonnage regulier :ou 
non peut s'analyser comme la somme lineaire de plusieurs phenomeries 
independants, le variogramme du phenomene global correspondant a la 
somme lineaire des variogrammes de chacun des phenomeries 
20 independants qui le composent. 

Classiquement, on decompose le variogramme correspondant aux 
donnees experimentales mesurees en une somme de variogrammes 
modelis6s et on deduit des donnees experimentales et des modelisations 
choisies pour les variogrammes eiementaires de la decomposition, les 
25 fonctions el6mentaires qui d§composent la fonction aleatoire correspondant 
au phenomene global. 

II est ainsi possible d'extraire d'une cartographie de donnees 
sismiques du type de celles illustrees sur la figure 1 (par exemple, des 

a 

donnees experimentales brutes) d'une part le bruit blanc present dans ces 
30 donnees (figure 2a), d'autre part un bruit correspondant a des raies lineaires 
(figure 2b) et enfin des donnees filtrees, debarrassees de Pun et Tautre de 
ces deux bruits (figure 2c). 
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Les calculs de determination par krigeage des valeurs des fonctions 
elementaires decomposant une fonction aleatoire globale sont en eux- 
memes classiquement connus par rhomme du metier. 

On pourra par exemple a cet egard se referer aux articles et 
5 publications mentionnees dans la bibliographie jointe a la fin de la presente 
description. 

De facon tres generate, on determine la valeur d'une fonction 
elementaire intervenant dans la decomposition de la fonction aleatoire 
globale comme etant une combinaison lineaire de valeurs experimentales 
10 pour des points d'un voisinage immediat du point considere, ces valeurs 
experimentales etant affectees de coefficients de ponderation. 

En d'autres termes, si Ton considere qu'une fonction Z(x) se 
decompose en la somme de fonction elemejitaire Y u (x), 

ce. qui se traduit par : 

15 



on estime la composante Y u (x) par 

20 r**(x)= ± x a z a 



ou a est un indice muet designant les points consideres autour du 
point x pour lequel on cherche a determiner la valeur estimee Y u *(x), Z x 
etant la valeur au point x, N etant le nombre de ces points. 
25 On d§montre que les coefficients de ponderation Xa verifient 

Tequation : 



C C i** 



^01 



c u 

\yoN j 
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ou Tindice 0 designe le point pour lequel on cherche a determiner 
une estimation, les valeurs C 0 " a C% N etant des valeurs de covariance 

calculees a partir du modele u correspondant a la composante Y u (valeurs 
de la fonction de covariance pour les distances entre chaque point de 
5 donnee et le point a estimer), les valeurs Cg etant des valeurs de covariance 
calcul6es en fonction de la mod6lisation choisie pour le variogramme de la 
fonction a estimer (valeurs de la fonction de covariance pour les distances 
entre les points de donnees). 

Ces coefficients de pond^ration X x u se determinent done par simple 

10 inversion des matrices de covariance. 

PROBLEME(S) POSES PAR L'ETAT DE LA TECHNIQUE - PRESENTATION DE 
UINVENTION '< . 

15 Une des difficultes des techniques d'analyse Krigeante connuesria 

ce jour tient en ce qu'elles necessitent d'utiliser des modelisations des 
fonctions de covariances. 

Uinteret d'utiliser ces modelisations est qu'elles permettent de 
disposer de matrices qui sont d6finies positives et inversibles. 

20 On comprend toutefois que ces techniques de filtrage, si elles 

donnent de bons rdsultats, sont fortement d^pendantes de Texpertise 
individuelle de la personne qui choisie les modelisations des differents 
variogrammes. 

II y a IS une source d'erreurs, qui empeche que ces techniques 
25 puissent etre utilisees par des personnes qui n'en sont pas des specialistes. 

En outre, le choix des modelisations est egalement a Torigine de 
pertes de temps importantes en production. 

Un but "de I'invention est de pallier cet inconvenient et de proposer 
une technique de filtrage par analyse krigeante qui puisse etre mise en 
30 ceuvre de fagon automatique ou quasi-automatique. 

Elle propose & cet effet un p Procede pour le filtrage d'au moins 
deux series de donnees representatives des valeurs d'au moins un 



m m 

4 

parametre physique sur une meme zone, caracterise en ce qu'on determine 
ie variogramme croise de ces series de donnees et on resout I'equation de 
co-krigeage qui en resulte pour en deduire de facon automatique une 
estimation de la composante commune a ces series de donnees. 

Elle propose egalement un procede pour le filtrage d'au moins une 
serie de donnees representative des valeurs d'au moins un parametre 
physique sur au moins une zone, caracterise en ce qu'on identifie une 
modelisation d'une composante de variabilite spatiale de son variogramme, 
on soustrait cette modelisation du variogramme experimental et on resout 
I'equation de krigeage correspondant a ce variogramme differentiel pour en 
deduire une estimation de la composante de variabilite correspondante sur 
la serie de donnees. 

PRESENTATION DES FIGURES 

D'autres caracteristiques et avantages de I'invention ressortiront 
encore de la description qui suit, laquelle est purement illustrative et non 
limitative et doit etre lue en regard des dessins annexes sur lesquels : 

- les figures 1 et 2a, 2b et 2c, deja analysees, illustrent un exemple de 
cartographie de donnees sismiques et de decomposition par analyse 
krigeante correspondante ; 

- les 3a, 3b illustrent deux cartographies d'une meme zone, obtenues a 
partir d'acquisitions reaiisees en deux temps differents ; 

- la figure 4 est une cartographie de la composante commune aux 
cartographie des figures 3a et 3b ; 

- les figures 5a et 5b et les figures 6a et 6b sont des cartographies des 
composantes autres que la composante commune des cartographies 
des figures 3a et 3b ; 

- les figures 7a a 7c sont des graphes illustrant la repartition des erreurs 
dans le cas respectivement d'un filtrage standard, d'un filtrage par 
krigeage factoriel classique et d'un filtrage conforme a celui propose par 
I'invention (co-krigeage factoriel). 
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DESCRIPTION D'UN OU PLUSIEURS EXEMPLES DE MISE EN a=l/Vf?r= 
DE L'INVENT/ON 

Filtraqe automatiaue 

On se place ici dans le cas ou on dispose de deux cartographies 
obtenues pour une meme zone avec des donnees acquises par exempre a 
des instants differents ou encore pour des attributs sismiques differents. 

Ces deux cartographies sont par exemple du type de celle illustree 
sur les figures 3a et 3b. 

Dans ia suite du texte, on designera par 21 et Z2 les deux fonctions 
correspondents a ces deux series de donnees. 

II est propose de decomposer chacune de ces deux fonctions en ia 
somme de leur composante commune et de residus orthogonaux. s 

A cet effet, on determine a partir des deux series de donnees dont 
on dispose un variogramme croise, dont les valeurs sont : 

Y12 (h)= -1 ]T (Z1 (x) - Z1 (x+h))(Z2(x) - Z2(x+h)) 

ou x et x+h designent les couples de points pris en consideration selon la 
direction et pour la distance h pour lesquels on detennine la valeur du 
variogramme et 

oil N est le nombre de ces couples de points pour cette direction et cette 
distance. 

Connaissant ce variogramme croise, on determine ensuite une 
estimation de la fonction qui lui correspond, laquelle verifie : 

Z*i 2 <x) = £ ^Zl+2 X}ZJ 

oil les coefficients de ponderation \* a et X 2 p se determinent par inversion de 
I'equation de co-krigeage : 
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On notera que la matrice qui apparaTt dans cette Equation a 
I'avantage d'etre inversible sous certaines conditions de calcul. 
5 De cette fagon, en utilisant les covariances exp6rimentales, on 

decompose de fagon automatique les deux variables correspondant aux 
deux series de donnees initiates en une composante commune et deux 
composantes residuelles orthogonales. La regularity des donnees fait que 
la covariance experimental est connue pour toutes les distances utilisees, 
10 aucune interpolation n'est necessaire, done la matrice est definie positive >; 

La fonction alors obtenue est une estimation de la composante 
commune aux deux series de donnees. 

La figure 4 represente une estimation de composante commune 
obtenue a partir des donnees correspondants aux cartographies des figures 
15 3aet3b. 

On comprend qu'il est particulierement avantageux dans de 
nombreuses applications et notamment en sismique 4D de disposer de 
cette composante commune. 

Elle permet notamment d'6liminer les artefacts d'acquisition non 
20 repetables 

- a partir de grilles d'attributs sismiques et notamment par exemple a 
partir d'amplitudes RMS dans un intervalle 

- a partir de grille de temps sismiques et par exemple & partir de temps 
d'6venements sismiques 

25 - a partir de volumes de vitesses sismiques et par exemple a partir de 
pointes automatiques de vitesses. 

Elle peut egalement etre utilisee en sismique simple pour eliminer 
des artefacts d'acquisition non repetables 
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- notamment a partir de grilles d'attributs sismiques calcules sur les 
classes d'incidence consecutives ou encore 

- a partir de grilles d'attributs sismiques calculus sur des volumes issus de 
sommations partielles ou d'ondes converties. 

5 Par ailleurs, une fois cette estimation de la composante commune 

d6termin6e, il est possible de determiner les composantes residuelles 
correspondant a la difference entre les donnees initiates et cette 
composante commune estim§e. 

Ces composantes residuelles peuvent elles-memes etre 
1 0 decomposees par analyse krigeante. 

C f est ce qu'illustrent les figures 5a, 5b et 6a, 6b, qui sont des 
cartographies des bruits blancs et des bruits de rates lineaires estimes de 
cette fagon pour Tune et I'autre des deux series de mesures representees 
sur les figures 3a et 3b. 
15 Des exemples de mesures d'erreurs obtenues par filtrage standard, 

par filtrage par analyse krigeante puis par filtrage par analyse cokrigegnte 
(ou krigeage multivariate) sont illustres sur les graphes des figures 7a a ; 7c. 

On comprend a la lecture de ces figures que le filtrage par analyse 
cokrigeante permet des dispersions bien moindres que dans les cas des 
20 filtrages classiques ou de filtrages par analyse krigeante et donne de 
meilleurs resultats. 

On a decrit dans ce qui precede un exemple mis en oeuvre avec 
deux series de donnees. 

Comme on le comprendra aisement, le proced6 propose peut 
25 egalement etre mis en oeuvre de la meme fagon avec un nombres de series 
de donnees (campagnes) plus important. 

Filtrage semi-automatique 

30 Ce deuxieme mode de mise en oeuvre permet egalement une 

decomposition simpliftee dans le cas ou Ton ne dispose que d'une seule 
serie de mesures (fonction S1). 
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II suppose que Ton dispose au prealable de la modelisation d'une 
composante 9m du variogramme experimental 0. 

Cette modelisation que Ton connaTt au prealable est par exemple la 
modelisation d'une composante independante de la geologie du sous-sol : 
5 bruit blanc, stries, etc. 

Connaissant cette modelisation de la composante 9m, on en deduit 
le variogramme residuel correspondant a la difference entre le variogramme 
experimental et cette composante Gm. 

On met alors en oeuvre une analyse krigeante pour determiner d'une 
10 part la composante modelis6e Sm et d'autre part, a partir du variogramme 
residuel, le residu orthogonal R1 tel que : 

S1 = Sm + R1 

Cette decomposition automatique permet de filtrer les anomalies 
d'acquisition qui presentent une coherence spatiale aisement identifiable et 
15 modelisable, comme les stries paralleles aux cables observees sur les 
amplitudes et les temps en sismique marine. 

Les exemples donnes ci-dessous ont ete d£crits dans le cas de 
donn6es sismiques. 

On comprend toutefois aisement que llnvention peut s'appliquer de 
20 fa<?on generale a tout autre type de donnees, notamment a tout autre type 
de donnees cartographiques 
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REVENPICATIONS. 

1. Procede pour le filtrage d'au moins deux series de donnees 
representatives des valeurs d'au moins un parametre physique sur une 

5 meme zone, caracterise en ce qu'on determine le variogramme croise de 
ces series de donn6es et on r6sout Pequation de co-krigeage qui en resulte 
pour en d6duire de fagon automatique une estimation de la composante 
commune a ces series de donnees. 

2. Procede selon la revendication 1, caracterise en ce que Ton 
10 soustrait des donnees initiates Testimation de la composante commune 

ainsi determine, afin de deduire les residus orthogonaux pour les 
differentes series de donnees initiates. 

3. Procede selon Tune des revendjcations precedentes, caracterise 
en ce qu'on met en oeuvre une analyse krigeante pour decomposer lesdits 

15 residus orthogonaux. 

4. Procede pour le filtrage d'au moins une serie de donnees 
representative des valeurs d'au moins un parametre physique sur au moins 
une zone, caracterise en ce qu'on identifie une modelisatioa d'une 
composante de variabilite spatiale de son variogramme, on soustrait cette 

20 modelisation du variogramme experimental et on resout I'equation de 
krigeage correspondant a ce variogramme differentiel pour en d6duire une 
estimation de la composante de variabilite correspondante sur la serie de 
donn6es. 

5. Proc6de selon Tune des revendications precedentes, caracterise 
25 en ce que la s£rie de donnfees ou au moins Tune des series de donnees est 

representative de donnees sismiques. 
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